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The origin of the breaking of conventional linear fcj_-factorization for hard processes in a nuclear 
environment is by now well established. The realization of the nonlinear nuclear fcj_-factorization 
which emerges instead was found to change from one jet observable to another. Here we demonstrate 
how the pattern of nonlinear fcj_-factorization, and especially the role of diffractive interactions, in the 
production of dijets off nuclei depends on the color properties of the underlying pQCD subprocess. 
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PACS numbers: 13.97.-a, 11.80La,12.38.Bx, 13.85.-t 

The fundamental point about the familiar perturbative 
QCD (pQCD) factorization theorems is that the hard 
scattering cross sections are linear functionals (convolu- 
tions) of the appropriate parton densities in the projectile 
and target. A consistent analysis of forward hard dijet 
production in deep inelastic scattering (DIS) off nuclei 
revealed a striking breaking of fc^-factorization con- 
firmed later on in the related analysis of single-jet spectra 
in hadron-nucleus collisions Q. Namely, following the 
pQCD treatment of diffractive dijet production 0,0] one 
can define the collective nuclear unintegrated gluon den- 
sity such that the familiar linear k± -factorization (see e.g. 
the recent review ) would hold for the nuclear structure 
function F 2 a(x, Q 2 ) and the forward single-quark spec- 
trum in DIS off nuclei because of its special abelian fea- 
tures. However, the dijet spectra and single-jet spectra 
in hadron-nucleus collisions prove to be highly nonlinear 
functionals of the collective nuclear gluon density. Fur- 
thermore, the pattern of nonlinearity for single-jet spec- 
tra was shown to depend strongly on the relevant par- 
tonic subprocess 2] . Our conclusions on the breaking of 
linear fcj_-factorization for hard scattering off nuclei have 
recently been corroborated by other authors [a, Q- In 
this communication we investigate the nonuniversality 
aspects of nonlinear nuclear fc^-factorization for differ- 
ent dijet excitation processes. We also comment on the 
process-dependence of the significance of diffractive final 
states. Our results shed certain light on to which extent 
hard processes in a nuclear environment can be described 
entirely in terms of the classical gluon field of a nucleus 

Here we extend the analysis |9(] of the excita- 
tion of heavy flavor and leading quark jets in DIS, 
7* — > QQ, to the excitation of open charm (or hard 
quark-antiquark dijet) and gluon jets in subprocesses 
9*9 ~ * QQi 1* 9 19 which are of direct relevance 
to the the large (pseudo) rapidity region of proton-proton 
and proton- nucleus collisions at RHIC. Our treatment 
is applicable when the beam and final state partons 
interact coherently over the whole nucleus, which at 
RHIC amounts to the proton fragmentation region of 
x = M 2 jj/2m p E b <x A = l/R A m p « 0.1A -1 / 3 , where 
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FIG. 1: The S-matrix structure of the two-body density ma- 
trix for excitation a — > be. 



Ra is the radius of the target nucleus of mass number 
A, E a is energy of the beam parton a in the target rest 
frame and m p is the proton mass [Tot ITU . 

To the lowest order in pQCD, all the above processes 
are of the general form ag — > be and, from the labora- 
tory frame standpoint, can be viewed as an excitation of 
the perturbative 1 6c) Fock state of the physical projec- 
tile | a) by one-gluon exchange with the target nucleon. 
In the case of a nuclear target one has to deal with the 
non-abelian intranuclear evolution due to multiple gluon 
exchanges which are enhanced by a large thickness of 
the target nucleus. The derivation of the master formula 
for di jet sp ectrum, based on the technique developed in 
HI llil Il3l| . is found in Q ; here we only reproduce the 
main result: 



dcr{a* — > be) 
dzbd 2 p b d 2 p c 
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d 2 b b d 2 b c d 2 b' b d 2 b' c 



(2tt)4 

x exp[-ip b {b b - b' b ) - ip c (b c - b' c ) 
xV(z b ,b b -b c )y*(z b ,b' b -b' c ) 



x{S^ b (b' b ,b' c ,b b ,b c ) + Sk 2 J(b',b) 
-S^(b,b' b X)-S^ c (b',b b ,b c )}.(l) 

If b a = b is the impact parameter of the projectile a, 
then b b = b + z c r, b c = b — z b r , where z b>c stand for 
the fraction the lightcone momentum of the projectile a 
carried by the partons b and c. All S^ n ' describe a scat- 
tering on a target of color-singlet systems of partons, as 
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indicated in Fig. ^ and all the dijet spectra are infrared 
finite observables. The and are readily calcu- 
lated in terms of the 2-parton and 3-parton dipole cross 
sections 0, Q, ^| , general rules for the multiple scat- 
tering theory calculation of the coupled-channel S 1 - 4 ' 1 are 
found in and need not be repeated here. 

The fundamental quantity of the color dipole formal- 
ism is the qq dipole cross section 0, 0] 



<Jqq{x.r) = 
f(x,K) = 



4iras(r) 



(2) 



where T(x,k 2 ) = dG(x, n 2 )/d\ogK 2 is the unintegrated 
gluon density in the target nucleoli. It furnishes a 
universal description of the proton structure function 
F2p(x, Q 2 ) and of the final states. For instance, the linear 
fc^-factorization for forward dijet cross section reads (for 
applications, see and references therein) 



2(2ir) 2 d<r N ( 1 * -> QQ) = 
dzd 2 p_d 2 A 

/(as, A) \*(z,p_) - *(z,p_ - A) | 



(3) 



where ty(z,p) is the qq wave function of the photon and 
A = p + +p_ is the jet-jet decorrelation momentum, and 
p , z = Z- refer to the Q— jet. 

The principal issue is the pQCD description of the dijet 
cross section for nuclear targets. First place, one needs 
the collective nuclear gluon density (j)(b,XA,K) per unit 
area in the impactparameter plane. As a starting point, 
one can define it [l|, |^ in terms of the qq nuclear profile 
function: 



r 2 A[b,(Tqq(x,r)} = 1-exp 



1 



a(x,r)T(b) 



d 2 Ktp(b,x,n) 1 — exp(inr) (4) 



It satisfies the sum rule / d 2 K(j}(b,x,K) = 1 — S a bs(b) , 
where S a b s (b) = exp[— ^ao(x)T(b)] and <Jo( x ) = 
J d 2 nf(x,n) is the dipole cross section for large dipoles 
and T(b) is the optical thickness of a nucleus. The 
explicit expansion for <j>(b, xa,k) in terms of the col- 
lective gluon density for j overlapping nucleons in the 
Lorentz-contracted nucleus [1(1 ] , its saturation properties 
at small n and the saturation scale QA(b, x) are found in 
P, 0] . Here we only emphasize that the inclusive spec- 
trum of leading quarks from the excitation 7* — > QQ 
off nuclei satisfies the same linear /cj_-factorization in 
terms of <p(b, xa, n) as its counterpart for free-nucleon 
target in terms of f(x,n), see Eq. ©. I.e., all Initial 
and Final State distortions of the spectrum of leading 
quarks are reabsorbed into the collective nuclear gluon 
density; such an abelianization only holds for the color 



singlet projectile. We shall also make use of <I>(b, x, k) — 
, xa, k). 

The i-channel pQCD gluon exchange leaves the target 
nucleon debris in the color excited state. In the case of 
nuclear targets one must distinguish the coherent diffrac- 
tive and truly inelastic processes. In the former the be 
color dipole is in the same color state as the projectile 
a and coherent diffraction aA — > (be) A with retention 
of the target nucleus in the ground state is possible. It 
gives rise to exactly back-to-back dijets, i.e., the diffrac- 
tive contribution is oc (A) (for finite -size nuclei the 
A-dependence is controlled by a slightly modified nuclear 
form factor with the width A 2 < R 2 A , see Ref. 0). 

At large- N c considered here, excitation from the color- 
singlet to color-octet dipoles in truly inelastic DIS is sup- 
pressed oc l/iV c , see the matrix element a is in Ref. 
but this smallness is compensated for by a large num- 
ber of the octet dipole states (at arbitrary N c one must 
speak of the adjoint and fundamental representations, re- 
ferring to them as octet and triplet states must not cause 
any confusion). That is a reason behind the unitarity- 
induced coherent diffraction making ~ 50% of the total 
DIS off heavy absorbing nuclei 01 • I n DIS off nuclei, the 
dipoles first propagate as color-singlets, then at depth (3 
in a nucleus excite into the octet state and further color 
exchanges at the remaining depth [/3, 1] only rotate the 
dipole within the the octet state. With inclusion of the 
diffractive component the dijet spectrum from excitation 
7* — > QQ in DIS off nuclei equals L] 

(2v) 2 do A (T ^ QQ) i m , f 1 AR f , 2 ,2 

d 2 bdzd 2 p_d 2 A = 2 T{b) J d(3 J d Kld 
xf(x A , - (3, b, x A , A - Ki - k)$(1 - /3, b, x A ,Ki) 



*(/3; z,p - ki) - *(/?; z,p - m - k) 



+5^(A)\^(l;z gi p_)-^(z a ,p_)\ , (5) 

where 3>(/3, b, x, k) is the collective nuclear glue for the 
slice (3 of a nucleus defined by 



exp 
and 



±(3a(x,r)T(b) 



Jd 2 . m b,x,K)e M ^r) 



p_) 



J d 2 n$(l3,b 7 XA,K)^(z,p_ + k) 



is the wave function of the incident color-singlet dipole 
distorted by the coherent Initial State Interaction (ISI) 
in the slice (3 of a nucleus. The diffractive component 
is a quadratic functional of the collective nuclear glue. 
The first component in JSJ describes truly inelastic DIS. 
Here the slice (1 — (3) in which the dipole is in the color- 
octet state gives the Final State Interactions (FSI). The 
singlet-to-octet transition is described by the free-nucleon 
gluon density /(xa,h). In contrast to the free-nucleon 
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result ® the nuclear dijet spectrum is of fifth order in 
gluon field densities: a quartic functional of the collective 
nuclear glue for two slices of a nucleus and a linear one 
of the free-nucleon glue; it can not be described by the 
classical gluon field of the whole nucleus. In DIS the 
FSI looks as an independent broadening of the quark and 
antiquark jets. The nonlinear k±- factorization result (J5J 
must be contrasted to the free-nucleon spectrum JHJ; it 
entails nuclear enhancement of the decorrelation of dijets 
from truly inelastic DIS, the semihard dijets, \p±\ 2 < 
Q 2 A (b,XA), are completely decorrelated. 

The large- N c properties of excitation q* — > qg are sim- 
ilar to those of excitation 7* — > qq in DIS. The free- 
nucleon dijet cross section equals 

2(2%) 2 da N (q* ^ gq) _ 



dz g d 2 p g d 2 A 



^{z g ,P„ 



\V(z g ,p q - A) - 4<(z 9 ,p„ - z g A)\ 



(6) 



where now ^(z,p) stands for the wave function of the \qg) 
Fock state of the quark. Eq. © is simply the differential 
form of the single-jet spectrum derived in Ref. Q. 

The extension to nuclear targets is straightforward. 
The set of color singlet 4-parton states qgq'g' which en- 
ter the master formula includes |33), 66) and 1 1515) 
states (and their large- N c generalizations). The ampli- 
tude of excitation of the 66) and 1 1515) states from the 
initial state 1 33) is suppressed oc 1/N C , which is com- 
pensated for in the dijet cross section by the number of 
color states in 1 66) and |1515). At large N c one of the 
1 66) ± 1 1515) states decouples from the initial state 1 33) 
[l9j. The nuclear dijet spectrum takes the form 

(2n) 2 da A (q*^ qg) 1 f 1 ^ t ^^^^ 



d 2 bdzd 2 p g d 2 A 

x$(l - f3,b,XA, A - Ki - k)$(2 - (3,b,XA,K%) 

2 

x ^((3;z gi p g - Ki) - %(0;z g ,p g -m-K 
+<f>(b,x A , A) $(l;z g ,p g - A) - *(z„,p. 
+S^ 2 \A)S a bs(b) ^(l;zg,p g ) - ty(z g ,p g ) 



g,fg z g^-i 



(7) 



Here the third component in is the contribution 
from the coherent diffractive excitation of color-triplet 
qg dipoles, q* A — > (qg)A. It is suppressed by the nuclear 
attenuation factor which is a consequence of the initial 
parton q* being a colored one. The second term in Q 
can be associated with excitation of the color-triplet qg 
states. It looks like satisfying the linear fc^-factorization 
in terms of c/)(b, xa, A) but it does not: one of the wave 
functions, z g ,p g ), is coherently distorted over the 
whole thickness of the nucleus. Finally, the first compo- 
nent of the nuclear spectrum Q describes excitation of 



the color sextet and 15-plet qg states. The free-nucleon 
result JBJ is recovered to the impulse approximation. 

The latter contribution to the nuclear dijet spectrum 
is a fifth order functional of gluon densities and resem- 
bles strongly the truly inelastic dijet spectrum of (J3J) 
for DIS. As it was the case for DIS, the free-nucleon 
gluon density f{xA,i^) describes the excitation of the 
qg color dipole from the lower (triplet) to higher (sex- 
tet and 15-plet) color states. The ISI distortions are 
similar too, the principal difference is in the nuclear 
thickness dependence of the distortion factors in the 
second line of Eqs. (5) and (7): the asymmetric one, 
$(1-/3, b, xa,A-k 1 -k)^(2-(3, b,XA, Ki) for the frag- 
mentation of colored quark q* vs. the symmetric one, 
$(1 - /?, b, xa, A - ki - k)<I>(1 - /3, 6, xa, «i) in DIS. 
In DIS it describes equal distortion of the both outgo- 
ing parton waves by pure FSI, for the incident quarks q* 
in pA collisions it includes the ISI distortion of the in- 
coming wave of the colored quark q* . For the latter rea- 
son the input nonlinearity of this contribution is of still 
higher, sixth, order: the factor $(2 — (3, b, xa,k>i) which 
looks as if defined for the slice of nuclear matter of thick- 
ness (2 — f3) is, as a matter of fact, a convolution $(2 — 
/?, b, x A , A) = J d 2 K${l, b, x A , k)$(1 - A b,x A , A - «), 
which combines the effects of distortion of the incident 
quark wave over the whole thickness of the nucleus and 
of the produced quark in the slice (1 — ff) of the nucleus. 

The large- N c properties of heavy flavor excitation via 
g* — ► QQ are quite different. The free-nucleon dijet cross 
section is again the differential form of the single-jet spec- 
trum derived in 0]: 

2(2Tr) 2 da N (g* -> QQ) = 
dzd 2 p_d 2 A 

= f(x A , A) \\V{z,p_) - *{z,p_ - zA)\ 2 



•|*(«,p_ - A) -¥(*,p_ - 



;A)| 5 



(8) 



Because one starts with the color-octet QQ dipole, at 
large- N c the intranuclear interactions are color rotations 
in the space of octet QQ states. Transitions to the color- 
singlet cc dipoles are suppressed and the non-abelian evo- 
lution of the QQQ'Q' state becomes the single channel 
problem. The coherent diffraction excitation, in which 
the initial and final color states must be identical, is like- 
wise suppressed. The resulting nuclear dijet cross cross 
section equals 

{2ix) 2 da A {g* -> QQ) 



zA)\ 2 



dzd 2 p„d 2 bd 2 A = J d *mb,x A ,K) 
x$(l;b,XA,A - K)\ i $>(z,p_ — k) — ^{z,p_ 
= S abs {b)cl)(b,XA,A) 

x j|\l/(z,p_) - *(z,p_ - -zA)| 2 + 
+ |*(z,p_ - A) - *(z,p_ - zA)| 2 | 
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+ / d 2 Kcj){b, xa, xa, A — k) 

x\^(z,p_ - k) - *(z,p_ - zA)| 2 . (9) 

This result for the dijet spectrum © is precisely the dif- 
ferential version of the single-quark spectrum, Eq. (31) 
of Ref . Q] , if in the nonlinear term one makes an identifi- 
cation A = m + K2- It satisfies the quadratic- nonlinear 
fc_i_-factorization in contrast to the quartic-nonlinear one 
for the leading quark-antiquark dijets in DIS and the qg 
dijets from q* — > qg. Although only the collective glue for 
the whole thickness of the nucleus enters the nonlinear 
term in (J3J, and kincmatically it looks as if corresponding 
to the subprocess g* g\AgiA QQ with two uncorrelated 
collective nuclear gluons gA, the emerging combination of 
wave functions can not readily be associated with specific 
Feynman diagrams in terms of collective nuclear gluons 
gA- Different collinear contributions corresponding to 
poles of the wave functions ^{p_ + Ki) can readily be 
identified ( 0, 13| , for the related discussion see also |2cj): 
such an analysis goes beyond the scope of the present 
communication and will be addressed elsewhere. 

Our findings can be summarized as follows: The non- 
linear Ax-factorization relations and © for the dijet 
spectrum with two extreme color excitation properties, 
are our main new results. Such explicit representation for 
the dijet spectra is not contained in previous works on the 
subject [flLl- These examples show clearly that the non- 
linear fcj_-factorization is a generic feature of the pQCD 
description of the dijet production in a nuclear environ- 
ment. We established how at large iV c the pattern of non- 
linearity depends on color properties of the relevant QCD 
subprocess: (A) excitation of dijets in higher color rep- 
resentations from partons in a lower representation typi- 
cally gives rise to the fifth, or even sixth, order nonlincar- 
ity of of the dijet spectrum in gluon fields, while (B) in 
the processes starting from already higher-representation 
partons inelastic interactions can be viewed as color ro- 
tations within the same representation and the nonlin- 
earity will be quadratic one. A part of the nonlinearity 
comes from the free-nucleon gluon density which emerges 
in all instances of excitation of higher color representa- 
tions (see also the related discussion of the l/(iV| — 1) 
expansion in Ref. Q). The coherent diffraction is not 
suppressed by large N c in the class-A reactions and is 
color-suppressed in the class-B reactions; this is also a 
new observation. Still another feature of the class-A re- 
actions is a contribution from dijets in the same color 
representation as the incident parton. Within our non- 
linear fcj_-factorization, all the dijet spectra are explicitly 
calculable in terms of the collective nuclear glue of Eq. 
(@J. However, in the class-A reactions this collective glue 
must be evaluated for different slices of a nucleus and 



enters ISI and FSI effects in quite a distinct way. The 
condition, x < xa ~ 0.1 • A" 1 / 3 , restricts the applicabil- 
ity domain of our formalism to the proton hemisphere of 
pA collisions at RHIC; the required coherency condition 
does not hold for the mid-rapidity dijets studied so far 
at RHIC 
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